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Abstract 

We derive robust model-independent bounds on Dark Matter (DM) annihilations 
and decays from the first year of Fermi 7-ray observations of the whole sky. These 
bounds only have a mild dependence on the DM density profile and allow the follow- 
ing DM interpretations of the PAMELA and Fermi e excesses: primary channels 
or e + e~e + e~ . An isothermal-like density profile is needed for 
annihilating DM. In all such cases, Fermi 7 spectra must contain a significant DM 
component, that may be probed in the future. 
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1 Introduction 



Recently the Fermi collaboration released the first sky map of 7-rays up to energies below a 
few hundred GeV [1, 2, 3]. From a particle physics point of view its main interest resides in 
the possible presence of a DM signal over the astrophysical background. Excitingly enough, 
although no clear excess is present, a few theorists claim possible hints [4]. In this paper we do 
not address these issues, that will need a full understanding of the data (released photon data 
still contain a non-negligible contamination from mis-identified hadrons at energies around and 
above 100 GeV), a proper modeling of the astrophysical backgrounds, subtraction of identified 
point-sources, and maybe more statistics. 

Here we take a different approach: using the available data we derive robust bounds on 
DM annihilations and decays, by only demanding the DM-induced 7-ray flux to be below 
the observed flux. Since DM is neutral, 7's are only produced at higher order in the QED 
coupling from various processes: we only include those contributions that can be computed in 
a model-independent way. 

We do not attempt to subtract the astrophysical and instrumental backgrounds, nor point 
sources: many of these subtractions can be performed only by assuming some astrophysical 
model, reducing the robustness of the bounds. Whenever possible, we will also comment on 
and try to quantify the residual uncertainties of our results. Progress on understanding the 
backgrounds can only render our bounds more stringent, presumably by order one factors. A 
Dark Matter signal could be lurking just below our bounds. 

An interesting application of our results is checking whether the DM interpretations of the 
e excesses observed by PAMELA [5], Fermi [6] and ATIC [7] give an Inverse Compton (IC) 
photon flux compatible with 7-ray observations. Indeed, according to DM interpretations, the 
e excess should be present everywhere in the DM halo (rather than only locally, as if the e 
excess is due to a nearby astrophysical source such as a pulsar), giving rise to an unavoidable 
associated 7-ray signal: produced by DM loose essentially all their energy by Compton up- 
scattering ambient light, giving rise to a photon flux at the level of Fermi sensitivity. Being a 
diffuse signal, one can consider regions of the sky that have smaller astrophysical uncertainties. 
By doing so, we find that many DM interpretations of the e excesses are already excluded by 
this bound (future improvements will not change this conclusion); moreover, in some cases the 
expected effect is at the level of the observed flux (future improvements should allow to test 
it). 

The plan of the paper is the following: in Section 2 we present technical details of our 
analysis; in Section 3 we present results for DM annihilations and in Section 4 for DM decays. 
Results are summarized in the Conclusions. 

2 Technical introduction 

We consider DM annihilations (parameterized by the DM DM cross section av) or decays 
(parameterized by the DM decay rate T = 1/r) into the following set of primary DM particles: 

2e, 2fi, 2r, 4e, 4/x, 4r, 2q, 2b, 2t, 2h, 2W (1) 
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where 2e stands for e + e~ and 4e stands for VV, where V is an hypotethical new light particle, 
whose mass we take to be m ~ 1 GeV, that decays into SM final states, as e.g. in the model 
of [8]. We denote as M the DM mass. For the Higgs boson, we assume a mass = 115 GeV. 
Concerning all the details not specified here, we follow [9]. 

2.1 Astrophysics 

We consider the following Milky Way DM tentative density profiles p(r) [10]: 

(1 + r /r 2 )/(l + r 2 /r 2 s ) isothermal, r s = 5kpc 

(r /r)(l + r /r s ) 2 /(l + r/r s ) 2 NFW, r s = 20 kpc (2) 
exp(— 2[(r/r s ) a — (r /r s ) a ]/a) Einasto, r s = 20 kpc, a = 0.17, 

keeping fixed the local DM density p(r = r 8.5 kpc) = p = 0.3 GeV/ cm 3 . NFW and 
Einasto profiles are favored by iV-body simulations, isothermal-like profiles by observations of 
spiral galaxies [11]. 

When assessing the residual uncertainties of our bounds, we will also consider the possibility 
of a disk-like component for the DM [12]. For this "Dark Disk" we will assume a profile 
Pdisk(r,z) oc exp(r /r^) smh(z/zd) 2 and we will vary r<i between 5 and 15 kpc, Zd between 0.2 
and 10 kpc, and the fraction of DM in the Dark Disk at solar position from to 50%. 

Regarding the diffusion of in the Milky Way, we consider the min, med, max propagation 
models of [13] characterized by the following astrophysical parameters: 

Model 5 Kq in kpc 2 /Myr L in kpc 



min 


0.85 


0.006 


1 


med 


0.70 


0.011 


4 


max 


0.46 


0.076 


15 



(3) 



The diffusion coefficient K = KqE s is assumed to be constant inside a cylinder with height 2L 
centered on the galactic plane and radius 20 kpc, and infinitely large outside. 

When assessing the dependence of our bounds on the size of the diffusion zone, we will 
consider also the more realistic case in which the diffusion coefficient depends on the distance 
from the galactic plane [14], K(E,z) = KqE s exp(|z|/ Zh). where Zh effectively determines the 
thickness of the diffusion zone, that dies off gradually instead of turning off abruptly at \z\ = L. 

2.2 Fitting the Fermi sky map 

We divide the Fermi 7-ray sky, parameterized by galactic longitude £ and latitude b (£ = b = 
corresponds to the Galactic Center, GC) into several regions, depicted in Fig. 1. We extract the 
photon spectrum within each region from the first public release of the Fermi 7-ray data [1]. 
The details of the extraction are summarized in the Appendix. The Fermi collaboration 
published the energy spectra below 100 GeV in a few regions [2] ; we checked that our procedure 
reproduces the Fermi results at low and intermediate energies. We do not subtract point 
sources, but we exclude the region most polluted by astrophysical sources, the galactic plane, 
by restricting ourselves to |6| > 5°. Since the signal we are seeking to bound is not uniformly 



P(r) 
Pe 
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Figure 1: Subdivision of the sky and example of a bound from each different region. In paren- 
thesis, the same bound is computed neglecting both the diffusion of e ± and the finite volume of 
the Milky Way diffusion halo. 



distributed on the sky, we consider a finer grid around the GC, where DM (but also astrophysical 
sources) concentrate. Conversely we use regions with increasingly larger areas towards the 
two poles, where the signal is expected to be smaller. Furthermore, we separately consider 
north-west, south-west, north-east and south-east regions, because the one less polluted by 
astrophysical sources will offer the best sensitivity. 

As the DM density profile is unknown, we do not know which region is most sensitive to DM. 
A robust constraint is obtained by demanding that the minimal computable DM-7 spectrum 
in each region, $p M j n energy bin i, does not exceed the measured 7 flux $° xp at 3a, for any 
energy bin and any region. Fig. 1 shows an example of the bounds on the DM annihilation 
cross sections for each region. Fig. 2a compares the Fermi data in the single region that gives 
the stronger bound for this particular model with the model prediction at its bets-fit point for 
the PAMELA and Fermi e ± excesses. 

Various regions give comparable bounds. Thereby one can do slightly better, still maintain- 
ing the absolute robustness of the bounds, by combining all regions in a global fit. We impose 
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Figure 2: Left: Fermi data compared with an example of best-fit DM annihilation signal. 
Photons above ~ 100 GeV can be still contaminated by hadrons. The dotted line shows Inverse 
Compton computed neglecting e diffusion and the finite volume of the diffusion halo. Right: 
Fermi preliminary extra- galactic data [3] compared with an example of best-fit DM decay signal. 



the 3<T bound, x 2 < 9, where 

x2 = ^ (ra^'-^W - - *n + £ , (4) 

where Q(x) = 1 if x > and otherwise; the sum runs over all angular and energy bins (with 
mean energy Ej)\ the \ 2 must be marginalized (minimized in Gaussian approximation) over 
the energy-scale free parameter e; the last term in the x 2 accounts for the 6e ~ 15% Fermi 
uncertainty on the energy scale; that only has a minor effect on the bounds. 

Typically, such global bound is a factor of few stronger than the bound obtained demanding 
that no single point is exceeded at more than 3a. Furthermore, there are conceptual advantages. 
The 'single point' bound depends on how we choose the energy and angular binning and can 
be fully dominated by the single bin where a downward statistical fluctuation happened in the 
total rate. On the other hand, the global fit does not depend on the binning (in the limit where 
it is dense enough); e.g. one can even split one bin into two coincident bins without affecting 
the global fit. 

2.3 Computing 7's from DM 

Since DM has no electric charge, 7-ray production from DM annihilations or decays occurs at 
higher order in the electromagnetic coupling from many different processes with comparable 
rates: i) bremsstrahlung from charged particles and 7r° decays; ii) virtual emission, iii) loop 
effects; iv) astrophysical processes involving other particles produced by DM. We only consider 
the following two sources of 7-rays that can be robustly computed in a model-independent way: 
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1. FSR7, i.e. Final State Radiation emitted by the primary DM annihilation or in subsequent 
decays. With a slight abuse of terminology, we will include in this contribution also the 
photons from hadronic decays. This gives photons with the largest E 7 ~ M. 

Primary channels such as 2r, 2q, 2W give rise to tt°s that decay as tt° — > 27 giving a 7 yield 
larger than channels (such as 2e or 2/x) that only produce 7 from bremsstrahlung. Models 
involving new neutral light particles give the smallest 7 yield [15, 16]. We do not consider 
electroweak bremsstrahlung, that cannot be computed in a model-independent way: one needs 
to know in which electroweak multiplets the DM lies. The large corrections found in [17] if 
M 3> Airv are only present for those decay or annihilation channels that arise thanks to a 
non- vanishing Higgs vev v. All the channels we consider can be realized, in appropriate DM 
models, as effective operators that do not involve the Higgs (e.g. as effective operators where a 
DM pair couples to a vector leptonic current, or to the W, Z field-strength squared), such that 
electroweak bremsstrahlung remains one higher order effect, as in [18]. 

At lower photon energies, the DM7 flux is dominated by the second source of 7-rays that 
we consider [9, 19] 

2. IC7, i.e. photons from Inverse Compon. DM gives rise to e ± that loose most of their 
energy by up-scattering galactic ambient light (CMB and starlight, partially rescattered 
by dust): this Inverse Compton e ± 7 — > e l±r j' process gives rise to 7's with energy Ey ~ 
E 7 (E e /m e ) 2 ~ 30 GeV. 

The IC energy loss process competes with energy losses due to synchrotron radiation in the 
galactic magnetic fields. The rates of these two processes are respectively proportional to w 7 (af) 
and to ub{x) = B 2 /2, the energy densities in photons and in magnetic fields. We take the 
values from [20, 21], and in particular we assume 

B(r,z) « ll//G-exp(-r/10kpc- |*|/2kpc). (5) 

In this case w 7 ^> ub everywhere: Inverse Compton is dominant, and it can be reliably computed 
as essentially all the e ± energy goes into IC, irrespectively of the precise galactic maps of m 7 
and of ub- 

Only large deviations from the maps we adopted can affect our bounds, so we comment 
about this possibility. Both galactic radiation and magnetic fields are better known for our 
neighborhood than for other regions of the Galaxy, with the magnetic field being the most 
uncertain. The most realistic worry is that magnetic fields in the Inner Galaxy might be intense 
enough that one there has ub ~ w 7 , weakening our bounds. We checked the dependence on 
the magnetic field uncertainties by varying the scales in Eq. (5) on which the galactic magnetic 
field changes both radially and vertically, while keeping its value at solar position fixed. We 
find that a factor of 2 variation in these scales changes the IC7 fluxes at high latitudes by less 
than 10 -T- 20% and in regions closer to the GC (|6| < 20°) by 60%. Since these regions are 
relevant in our global fit, our bounds can be relaxed at most by a factor of ~ 1.5 4- 2. 

In order to compute IC7, following [9] we take into account the diffusion of e ± , with charac- 
teristic diffusion length of about A ~ 1 kpc, such that photons from IC are smeared in a region 
with angular size A/r Q ~ 10°: even if DM annihilations are concentrated close to the GC, IC7 
are not. Thereby Fermi 7-ray data are more relevant than HESS data both for the energy 
range and for the angular range they observe. 
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Figure 3: FERMI full-sky bounds on Final State Radiation 'j-rays, for the DM annihilation 
modes indicated along the curves. 




Figure 4: FERMI full-sky bounds on Inverse Compton 'y-rays, for the DM annihilation modes 
indicated along the curves. 
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Figure 5: Example of how the Fermi IChf bound on o~v changes as function of the height L of 
the diffusion volume for different DM profiles. We here assumed DM annihilations into /x + /i~ 
with M = 1.3 TeV, but this plot would be almost the same for other DM models. 

3 DM annihilations 

Assuming DM annihilations, fig. 3 shows the Fermi all-sky global bounds on FSR7 as function 
of the DM mass M and the DM cross section av. Fig. 4 shows the corresponding bounds on 
IC7, assuming L = 4 kpc. 

In both cases the left (right) panels holds for the isothermal (NFW) DM profile. We see 
that the Fermi bounds only have a mild dependence on the DM profile: although we do not 
know where DM is, Fermi observed all the sky, so that it is no longer possible to hide DM 
with an appropriate density profile. 

The diffusion volume of is assumed to be a cylinder that extends away from the galactic 
plane up to \z\ < L. Fig. 5 shows a typical example of how the IC bounds depend on L and on 
the DM profile. The main result is that if L is as small as 1 kpc, DM can produce a significant 
fraction of its e outside of the diffusion volume, that negligibly contribute to the IC7 signal. 
Indeed they escape away, as the e mean free path is one or two orders of magnitude longer 
than the Milky Way, and the probability of entering into the diffusive halo is small. We also 
checked the importance of the diffusion zone thickness with the more realistic model described 
in Sect. 2.1, in which diffusion exponentially dies off on a length scale z^. Solving the diffusion 
equation in an ideally infinite volume, we find that the IC7 flux is closer to the large L « 15 
kpc case even for zu = 2 4 kpc. Therefore any realistic bounds should be much closer to the 
"max" case than to the ones with very small L, especially since the energy dependence of the 
diffusion coefficient plays only a sub-leading role in the IC predictions. 

The possibility of a thin diffusion volume will be relevant for our later discussion, so that it 
becomes important to settle this issue. As far as we know, it is disfavored by various arguments: 
a) global fits of charged CR propagation models favor L f» 4 kpc but values between 1 and 15 
kpc are considered [13]. b) abundances of CR with a life-time comparable to the diffusion 
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Figure 6: Bounds on DM annihilations into leptonic channels. The Fermi bounds are 
denoted as FSRy (continuous blue line) and IOy (red curves, for L = 1,2, 4kpc from upper to 
lower). Other bounds are described in the text; their labels appear along the corresponding lines 
only when these bounds are significant enough to appear within the plots. Cosmological freeze- 
out predicts av ~ 3 1CT 26 cm 3 / sec (lower horizontal band) and connections with the hierarchy 
problem suggest M ~ (10 -j- 1000) GeV. The region that can fit the e excesses survives only 
if DM annihilates into e's or fi 's and DM has an isothermal profile. All bounds are at 3cr; 
the green bands are favored by PAMELA (at 3cr for 1 dof) and the red ellipses by PAMELA, 
FERMI and HESS (at 3 and ha, 2 dof, as in [9]). 



time [22]; c) more realistic boundary conditions as described above; and presumably d) the fact 
that Fermi observes 100 GeV 7 rays also away from the GC suggests that L is not small. 

Fig.s 6 and 7 show again the Fermi bounds at 3a (the IC7 bounds is plotted for a few values 
of the height of the diffusion volume, L = 1, 2, 4kpc), together with the regions favored by the 
e ± excesses and with various other 3a bounds already considered in previous papers [15, 23, 9]: 

- The GC-7 (blue continuous curves) and GR-7 (dot-dashed blue curves) bounds refer to the 
HESS observations [24, 25] of the photon spectrum above ~ 200 GeV (so that it constrains 
FSR7 and heavier DM, rather than IC7 and lighter DM) in the 'Galactic Center' region 
(Vi 2 + b 2 < 0.1°) and in the 'Galactic Ridge' region (\£\ < 0.8° and \b\ < 0.3°). In these 
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Figure 7: Bounds on DM annihilations into non-leptonic channels. These channels can 
fit the PAMELA e + excess, but not the Fermi e + + e~ excess. The Fermi bounds on FSR'y 
exclude non-leptonic DM interpretations of the PAMELA e + excess, even for an isothermal 
DM profile. Non-leptonic branching ratios must be small. 



regions the DM density p(r) is uncertain by orders of magnitude, such that one gets strong 
bounds assuming NFW-like DM profiles and negligible bounds assuming isothermal-like 
profiles. 

- The dS-7 bound (dashed blue curves) refers to the HESS and VERITAS observations of 
various dwarf spheroidal galaxies [26, 27, 28]. 

- The v bounds (black curves) refer to the SuperKamiokande (SK) observations of neutrino 
from regions around the Galactic Center [29, 23, 9, 30]. 

- The GC-radio bound (red dashed curves) refers to radio observations of the Sgr A* black 
hole at the dynamical center of the galaxy and depends on the extremely uncertain local 
DM density [31, 15]. 
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- The cosmological CMB bound (red dashed curves) refers to the contribution St to the 
optical depth of CMB photons due to DM re-ionization of H and He 1 . 

In conclusion, DM interpretations of the e ± excesses survive only if DM annihilates into 2p, 
Ap or 4e and if DM has a quasi- const ant isotermal-like density profile. The GC-7 and GR-7 
bounds already disfavored solutions involving NFW or Einasto profiles [15, 9], but this needed 
extrapolating these profiles down to small scales not probed by iV-body simulations. Now an 
unseen excess would be present at larger scales where iV-body simulations are under control and 
favor these profiles. Furthermore, channels involving r are now disfavored even for an isothermal 
profile. In view of the FSR-7 FERMI bound, non-leptonic channels (fig. 7) can similarly at 
most have a small sub-dominant branching ratio, so that solutions involving Minimal Dark 
Matter or the supersymmetric wino [33] are now firmly excluded. 

The allowed solutions predict that a sizable fraction of the photons observed by Fermi 
around 100 GeV must be due to IC7 from DM e . The Fermi bound on IC7 becomes weaker 
if the diffusive volume of our galaxy is thin, L w lkpc (dotted red curves). 

Another possible way of weakening the bound is assuming that a fraction of the local DM 
density is stored in a dark disk component, such that FSR7 and (to a lesser extent) IC7 are 
moved towards the Galactic Plane, where the astrophysical 7 background is higher. However, 
to relax the conclusions on the DM profile, one needs this fraction to be large, of order unity, 
especially if the dark disk has a thickness z<i not much smaller than r . 

On the other hand, the Fermi bound can be made stronger subtracting from the 7 spectra 
the hadrons misidentified as 7 and the identified astrophysical point-like sources. 

4 DM decays 

Interpretations of the e ± excesses in term of DM decays (rather than annihilations) attracted 
interest because they were not in tension with 7-ray observations [35] (see also [36]). Indeed, 
the space-time density of DM decays is p/tM while the space-time density of DM annihilations 
is av(p/M) 2 /2. Thereby, HESS observations of 7-rays from the Galactic Center gave significant 
constraints on DM annihilations if r 3 p(r) 2 is large for r — » 0. On the contrary, even for a NFW 
profile, r 3 p(r) remains small such that DM decays were not significantly constrained. 

Fig. 8 shows, in the mass-lifetime plane, the new Fermi FSR-7 an d IC-7 bounds on DM 
decays, together with previous bounds from SK neutrino observations (for simplicity we do not 
plot the previous HESS 7 bounds, as they are now subdominant) and with the regions favored 
by interpretations of the e excesses in terms of DM decays. We see that such interpretations 
are now constrained. Also in the case of DM decays, the remaining viable channels are 2p, 
Ap or 4e. Channels like r's, producing tt° and other mesons decaying into photons, are now 

1 We adopt the computation by Cirelli et al. [32] (not performed for all DM channels we consider), who 
plotted the WMAP bound at la, St < 0.064. In our plots the region suggested by the PAMELA and Fermi e ± 
excesses is compatible with the CMB bound because we plot the WMAP bound at 3c, St < 0.094. The CMB 
bound does not depend on the local DM density pq, while all other curves actually constrain av p Q . Ref. [34] 
claims that pq is larger than the p© = 0.3 GeV/ cm 3 assumed here and close to 0.4 GeV/ cm 3 : in such a case 
the CMB bound would be relatively less stringent by a factor of 1.8. 
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Figure 8: Bounds on DM decays. In the upper rows we consider the leptonic channels that 
can fit the e excesses. In the lower row we consider the 'traditional' channels. 



12 



disfavored. Even in this case, a significant fraction of the photons observed by Fermi around 
100 GeV away from the Galactic Center should be due to DM. 

The bound denoted as 'exG-7' is obtained demanding that the cosmological 7 flux from 
DM decays does not exceed the extra-galactic isotropic flux observed by Fermi [3]. This 
cosmological flux is expected to be comparable to the galactic flux: 



^cosmo Pcosmo-^cosiU) 



^galactic 



(6) 



where p C osmo = ^dmPct ~ 1-3 10 6 GeV/cm 3 and -R C osmo ~ ^/H ps 13Gyr. 2 The isotropic 
cosmological 7 flux is 

d% _c_ r 1 e~ T a dN^E™ = Eja) 

dKf ~ An Jo a H(a) dVdtdE™ ^ ' 

where the first term inside the integral generalizes the usual line of sight integ rand ds, to the 
cosmological geometry described by the Hubble rate H(a) = H ^jQ\ + Q m /a 3 , as function of 
the scale factor a of the universe, which gives the E™ = E^/a redshift 3 . The last term is 
the usual space-time density of 7 sources, equal to T(p-DM{a)/M)dN y /dE^ in the case of DM 
decays. We can neglect absorption of 7, as the optical depth is r 1 below a few TeV. Eq. (7) 
can be directly applied to the computation of FSR7. 

The flux of photons generated by cosmological Inverse Compton scatterings of e from 
DM annihilations on Comic Microwave Background with energy density u 7 (a) = 7rT 4 /15 and 
spectrum dn^/dE = E 2 /n 2 j{e E l T — 1) at temperature T = T /a can be written as: 

rf$ 7 c 9cmiT Po f i da/ a 2 ff in 1 dn 7 dEf dE™ 

7 



dE, 32vr H M Jo ^n A + n m / a 3 JJ y e X dE™ E^ E 
where N e (E) = J™ dE' dN e /dE' and the function 



f lc = 2glng + (1 + 2q)(l - q) + ~ q) (9) 

z ± ~~\~ €q 

describes IC scattering / y(E^)e(E 1 / P-) — > e r y{E° ut = E 1 /a) at E™ ^> m e in terms of the dimen- 
sionless variables 

Ef mj T(l — e) 

E° ut lies in the range E^/E™ < e < T/(l + r). The non-relativistic (Thompson) limit corre- 
sponds to T < 1, so that e < 1 and < q < 1. 

Fig. 2b shows the Fermi isotropic 7 data (according to the preliminary analysis in [3]), 
compared to the cosmological flux generated by the best-fit DM — > t + t~ model. In such a case 



2 This bound was considered in [37, 38] in the case of DM decays and in [39] in the case of DM annihilations. 
In such a case the larger suppression, (/o cosmo / Po) 2 is counteracted by DM clumping in structures and galaxies, 
an effect which cannot be computed reliably. 

3 Thereby the cosmological signal allows to probe DM decay models (not considered in this paper) that 
give hard photon emission only at 400 GeV — 1 TeV, as these photons, above the Fermi reach, get partially 
red-shifted down to the Fermi range. 
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the main constraint comes from FSR rather than from IC. Other DM decay models that can 
fit the e ± excesses can have a much smaller FSR7, but have a very similar IC7. Thereby all 
such models predict an IC7 at the level of the Fermi observations. Depending on how Fermi 
extracted the isotropic component of their sky map, it might be correct to include in it not 
only the extra-galactic DM7 flux but also the galactic DM7 flux from the direction where it is 
minimal, thereby strengthening the bound [38] . 

5 Conclusions 

We presented robust model-independent bounds on DM annihilations and decays obtained 
demanding that the computable part of the DM-induced 7 flux be below the observed flux, 
as recently observed by Fermi in the full sky. Fig. 1 shows an example of bounds from the 
different regions of the sky we consider. Our results in fig.s 3 to 7 for DM annihilations and in 
fig. 8 for DM decays are based on a full-sky global fit, and thereby only have a mild dependence 
on the DM density profile. 

We show the bounds on Final State Radiation 7 separately from the bounds on Inverse 
Compton 7, as the latter can be weakened if one allows for significant variations in astrophysics 
with respect to the models we consider: i) increasing the magnetic fields until synchrotron 
energy losses dominate over IC; this is presumably allowed only in the inner regions of the 
Galaxy, and would reduce our full-sky IC bounds by up to 1.5 4- 2. ii) IC bounds can be 
reduced by a factor of few if the e diffusion zone is very thin (L ~ 1 kpc) and terminates 
abruptly, see fig. 5. iii) a Dark Disk component comprising ~ 50% of the local DM density can 
also be invoked to weaken these bounds. 

On the other hand, subtraction of astrophysical backgrounds (such as identified point-like 
sources) and of mis-identified hadrons, still present in the Fermi data we fitted, can only 
strengthen our bounds, presumably by a factor of few. 

Present data are enough to make progress on testing DM interpretations of the e excesses 
observed by PAMELA, Fermi, ATIC. At the light of the new Fermi 7 bounds, the remaining 
allowed DM interpretations involve DM annihilations or decays into fi + fi~ , VV — > fi + fi~ fi + fi~ 
or VV — > e + e~e + e~ primary channels, r's in the final staes are now disfavored even in the DM 
decay case. Moreover, for DM annihilation, a quasi-constant isothermal-like density profile is 
needed. This profile is not favored by DM simulations, that suggests that DM is concentrated 
around the Galactic Center. In such viable interpretation may be obtained e.g. 

assuming that DM annihilates into intermediate particles V with a lifetime longer than a few 
kpc [41] or that DM decays. 

Even in these cases, the expected DM signal is at the level of the observed flux so that it 
will be interesting to improve the sensitivity with forthcoming cleaner data and more statistics. 
According to [4], present Fermi 7 data already suggest the presence of a 'Fermi haze' excess 
with an angular and energy spectrum compatible with the expected IC DM excess. 
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A Extraction of the Fermi data 

We used the Fermi 7-ray data from [1]. We selected events from the 'diffuse' class, which 
have tighter cuts to reject the CR background. We considered the first 64 weeks of data (up 
to MET 280417908). After removing the Earth albedo, the data was binned in 0.5° x 0.5° in 
galactic latitude and longitude and in 16 logarithmically spaced bins in energy. We computed 
the effective area and exposure times and used them to convert the binned events into the 
photon flux. We then combined the bins in the larger areas of Fig. 1. We assigned systematic 
uncertainties according to [42] for the energy bins below 100 GeV: 10% at 100 MeV, 5% at 0.5 
GeV, progressively increasing to 20% at 10 GeV and above. 

In this analysis we considered also events around and above 100 GeV: while these are still 
significantly contaminated by CR and the systematics of the publicly available tools are not 
completely studied, we feel that they are still usable for setting bounds the way we proceed 
in this paper. Only the possibility that the currently available software may overestimate the 
instrument effective area at high energies would render too stringent our bounds. Although 
this effect would be (partially) balanced by the reduction in the flux once the data is further 
cleaned from CR contamination, we conservatively decided to increase to 50% the systematic 
uncertainty associated to all the data above 100 GeV. 
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